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Degeneracy of discrete energy levels of finite-length, metallic single-wall carbon nanotubes depends 
on type of nanotubes, boundary condition, length of nanotubes and spin-orbit interaction. Metal-1 
nanotnbes, in which two non-equivalent valleys in the Brillouin zone have different orbital angular 
momenta with respect to the tube axis, exhibits nearly fourfold degeneracy and small lift of the 
degeneracy by the spin-orbit interaction reflecting the decoupling of two valleys in the eigenfunctions. 
In metal-2 nanotubes, in which the two valleys have the same orbital angular momentum, vernier- 
scale-like spectra appear for boundaries of orthogonal-shaped edge or cap-termination reflecting 
the strong valley coupling and the asymmetric velocities of the Dirac states. Lift of the fourfold 
degeneracy by parity splitting overcomes the spin-orbit interaction in shorter nanotubes with a so- 
called minimal boundary. Slowly decaying evanescent modes appear in the energy gap induced by 
the curvature of nanotube surface. Effective one-dimensional model reveals the role of boundary on 
the valley coupling in the eigenfunctions. 

PACS numbers: 73.63.Fg, 73.22.Dj 


I. INTRODUCTION 

Metallic single-wall carbon nanotubes (m-SWNTs) are 
ideal one-dimensional (ID) conductors of nanometer to 
micrometer length. Due to the confinement in the finite- 
length, energy levels of electrons are quantized and the 
eigenfunctions show standing wave behaviorJ^ Fourfold 
degeneracy of the discrete energy levels observed in the 
tunneling spectroscopy measurements has been consid¬ 
ered as an intrinsic property of SWNTs reflecting the 
two non-equivalent, degenerate valleys at the K and K' 
points in the two-dimensional (2D) Brillouin zone (BZ) 
together with two spins degrees of freedomi^^— Recent 
measurements with ultraclean SWNTs have found fine 
structures of the order of sub-milli-electron-volt in tun¬ 
neling conductance spectra caused by spin-orbit inter- 
action^i2r2& which lifts the fourfold degeneracy by spin 
splitting in each valleyjii“~ On the other hand, some ex¬ 
periments show degeneracy behaviors other than above, 
such as gate-dependent oscillation of twofold and fourfold 
degeneracyj^^^— and many of them have been owed to 
extrinsic effects such as impurities. 

In our previous study^ we pointed out the asymmetric 
velocities in the same valley for the m-SWNTs because of 
the curvature of nanotube surface. That is, the left- and 
right-going waves in the same valley have different veloc¬ 
ities, > u^^(> 0), where is the velocity 

of left-going (right-going) wave at the K valley, and we 

have the relations ^ and v'^ ^ because 

of the time-reversal symmetry. At the same time, we 
claimed the formation of vernier-scale-like energy spec¬ 
trum, in which two sets of the energy levels with a con¬ 
stant energy separation between the levels have a simi¬ 
lar but not-exactly the same separation for each set, if 
the strong valley coupling occurs, in which each set of 
the wavefunction is formed from a left-going wave at one 
valley and a right-going wave at another valley. As the 
result of the quantization of the wavenumber in the axis 


direction, there are two different sets of equi-spaced dis¬ 
crete energy levels, hv^'^ t :/and /iu^V/Lnt, like 
the vernier scale showing two- and fourfold oscillations 
as observed in the experiments,— liSr— where Lnt is the 
nanotube length. On the other hand, for the case of 
valley decoupling, in which each wavefunction is formed 
from a left- and right-going waves in the same valley, 
the fourfold degeneracy and its lift by the spin-orbit in- 
teractionii^— would be observed. Thus, it is important 
to reveal the coupling of the two valleys as a function 
SWNT chirality or the boundary shape for understand¬ 
ing the degeneracy behavior. 

As is known that the particle-in-a-lD-box model can¬ 
not be directly applied to the m-SWNTs because there 
are two left-going waves and two right-going waves, in 
general, the ratios of these traveling waves in a standing 
wave are determined by microscopic conditions such as 
the chirality and the boundary condition. Previous cal¬ 
culations have shown the standing waves oscillating in 
the length scale of carbon-carbon bond for the armchair 
nanotubes, in which each standing wave is constructed in 
the condition of strong valley couplingOn the other 
hand for the zigzag nanotubes, the slowly oscillating 
standing waves for doubly degenerate levels can be con¬ 
structed in the condition of valley decoupling The 
SWNTs have been classified in terms of the boundary 
condition. Using generalized parameters, McCann and 
Fal’ko classified the boundary conditions for the Dirac 
electrons in the m-SWNTs By employing microscopic 
analysis on the boundary modes for the honeycomb lat¬ 
tice, Akhmerov and Beenakker showed that, except for 
the armchair edge, zigzag-type boundary condition^ in 
which two valleys decouple each other in the eigenfunc¬ 
tions, is applicable for general boundary orientation with 
a so-called minimal boundary, in which the edge has min¬ 
imum numbers of empty sites and dangling bonds, and 
these numbers are the samei^ Above theory^^S as well 
as assuming the slowly varying confinement potential,™ 
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supports the decoupling of the two valley as a typical 
case for the m-SWNTs except the armchair nanotubes. 

Here, we will show that the two valleys are strongly 
coupled in the chiral nanotubes classified into so-called 
metal-2 nanotubes (see 1 }TT)) with certain boundaries, 
as well as the armchair nanotubes. The effect of the 
strong coupling combined with the asymmetric veloci¬ 
ties appears as the vernier-scale-like spectrum. For the 
so-called metal -1 nanotubes;^ and for the metal -2 nan¬ 
otubes with the minimal boundary, it will be shown the 
fourfold degeneracy and its small lift by the spin-orbit 
interaction as the result of decoupling of two valleys. In 
addition, it will be shown that parity splitting of the 
valley degeneracy overcomes the spin-orbit splitting for 
shorter metal -2 nanotubes with the minimal boundary. 

In this paper, we mainly focus on the finite-length m- 
SWNTs in which the ends and the center have the same 
rotational symmetry. We will show that the degeneracy 
of discrete energy levels of m-SWNTs strongly depends 
on the chirality, boundary condition, length and the spin- 
orbit interaction. We will revisit and analyse the cutting 
lines with the point of view of the orbital angular mo¬ 
mentum, then will perform numerical calculation for an 
extended tight-binding model and analytical calculation 
for an effective ID model to investigate the degeneracy 
and the valley coupling. 

This paper is organized as follows. In m occurrence 
of the valley coupling is discussed by analyzing the cut¬ 
ting lines. In llllll numerically calculated energy levels 
by using an extended tight-binding model for metal -2 
nanotubes with a couple of boundaries are shown. In 
m an effective ID model for describing the valley cou¬ 
pling is derived and the microscopic mechanism of the 
valley coupling is analytically investigated for a couple of 
boundary conditions. Analytical forms of the discretized 
wavenumber are also given. The conclusion is given in 
IfVl In the Appendix detailed analysis on the long 
cutting line is given. In the Appendix [B1 numerical cal¬ 
culation for metal- 1 , armchair and capped metal -2 nan¬ 
otubes is given. In the Appendices [C] and |D1 we discuss 
detailed derivation and mode analysis of the ID model, 
respectively. In the Appendix [El relation on coefficients 
of standing waves between A- and B-sublattices under a 
boundary is given. 


II. CUTTING LINE 

Let us consider a SWNT defined by rolling up the 
graphene sheet in the direction of the chiral vector 
Ch = nai + ma 2 = {n,m), where n and m are inte¬ 
gers specifying the chirality of SWNT, ai and 02 are the 
unit vectors of graphene The m-SWNTs, which satisfy 
mod(2n -|- m, 3) = 0, are further classified into metal-1 
(dfl = d) or metal-2 {da = 3d), where d = gcd(n, 77 i) 
is the greatest common divisor (gcd) of n and to, dji = 
gcd(2n -|- TO, 2 to -|- n), and it has been known that the 
K and K' points sit on the center of the ID BZ for the 


metal-I nanotubes, while they sit on 1/6 and 5/6 posi¬ 
tions for the metal -2 nanotubesi^ 

In this section, we will show that the two valleys have 
different orbital angular momenta for the metal -1 nan¬ 
otubes, whereas they have the same orbital angular mo¬ 
mentum for the metal -2 nanotubes from analysis of the 
cutting line, ID BZ plotted in 2D fc-spacei^ Here the 
orbital angular momentum of the valley is that at the 
valley center [K (K') point], and is given by an integer 
specifying the cutting line passing through the K [K') 
point. The corresponding properties have been shown in 
the previous work numericallyi^ Here, we will show a 
proof of this property analytically. 

States on a cutting line represent ID wavevectors in 
the direction of nanotube axis with an orbital angular 
momentum with respect to the nanotube axis which cor¬ 
responds to a wavevector in the circumference direction 
of the SWNT. For a finite-length SWNT, the ID wavevec¬ 
tors are no longer good quantum numbers. If the bound¬ 
aries have the same Cd rotational symmetry around the 
nanotube axis with that of the SWNT, the orbital angu¬ 
lar momentum specified by a cutting line is a conserved 
quantity. In this case, an electron with a wavevector is 
scattered at the boundary to ID states with the same 
orbital angular momentum, that is, the scattering within 
the cutting line. 

Let us briefly review on the definition of the cutting 
lines for discussing their orbital angular momenta. There 
are arbitrary definitions for a complete set of the cutting 
lines as well as there are arbitrary definitions for 2D BZ 
of graphene as shown in Fig. |T] The detailed descrip¬ 
tion on the definitions of the cutting lines is found in the 
review article.— Instead of the conventional definition of 
the cutting lines with short segments the following def¬ 
inition of the cutting lines with long segments, 

with 

TT N TT N , , 

and /r = 0 ,--- ,d-I, ( 2 ) 

which is derived from the helical and rotational symme¬ 
tries^ is convenient to consider the properties under the 
Cd rotational symmetry. [See the long segments L 1 L 2 in 
Fig. [U defined by Eqs. © and ©■] Here the separation 
of cutting lines Ki = (— 12^1 + tib 2 )/N is perpendicular 
to the cutting lines and represents the discreteness of the 
wavevector in the circumference direction, and 


is the vector of short cutting lines in the conven¬ 
tional definition, 3^ where 61 = (27r/a)(I/\/3,1) and 
62 = (27r/a)(I/\/3,-I) are the reciprocal lattice vec¬ 
tors of graphene, ti and t 2 are integers defined by 
ti = (2to -I- n)/d}i, t 2 = —(2n + m)/dfi. T = 
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FIG. 1. Cutting lines for (a) {n,m) = (6,3), and (b) (7,4) 
SWNTs. In each figure, shadow areas show three different 
choices of the 2D BZ of graphene. S1S2 short segments denote 
conventional cutting lines, while L1L2 long segment is another 
definition of ID BZ. Here d = 3, drj = 3, T = \Ch\ ~ 

3y/7a and = 42 for (6,3) SWNT, and d = 1, d_R = 3, 
T = VsTa, \Ch\ = V^a and TV = 62 for (7,4) SWNT. 


a Y^3(r?'+^7?'+^w^/d/j is the ID nanotube lattice con¬ 
stant, N = 2{'n? -\- -I- nm)/dR is the number of A (B) 

atoms in the nanotube ID unit cell, a = 2.46A is the 
lattice constant of graphene. 

The rectangle defined by the two vectors dKi and 
NK 2 /d which surrounds the set of the longer cutting 
lines (see the vertically-long shadowed rectangle in Fig. 
[1]) is equivalent to the 2D BZ of graphene. The corre¬ 
sponding unit vectors in the real space are the vector 
Ch/d, and the vector given by i? = phO,i + qhO ,2 where 
Ph and qt satisfy the relation of mph — nqh = d^^^ 
The component of R in the direction of nanotube axis is 
expressed by 


for the given angular momentum p are represented in a 
single cutting line. Here, the orbital angular momentum 
of a state is defined by p in Eq. 0. 

It would be worthful to compare the present definition 
[Eq. (0)] with conventional definition of the cutting 
lines^ — tt/T < k < tt/T, and p = 0, - ■ ■ ,N — 1. [See 
the set of the short segments S 1 S 2 in Fig. m] In the 
conventional definition, every d cutting lines belongs to 
the same orbital angular momentum, that is, these cut¬ 
ting lines can mapped onto a single longer cutting line 
by translations with reciprocal lattice vectors. There¬ 
fore, the long cutting lines is convenient to consider the 
properties under the Cd rotational symmetry since an 
orbital angular momentum and a cutting line are one-to- 
one correspondence. 

Hereafter we focus on the m-SWNTs, in which there 
are cutting lines passing through the K and K' points. 
As proven in the Appendix [A] (and shown in Figs. [T] (a) 
and (b) as examples), the long cutting line L1L2 of p = 0 
passes through both K and K' points for the metal-2 
nanotubes, whereas cutting line passes only through K 
or K' points for the metal-1 nanotubes. The metal-2 
nanotubes are further classified into metal-2p and metal- 
2m by the conditions^^ 



for metal-2p, 
for metal-2m. 


(5) 


It is also shown in the Appendix [^ [and in Fig. [T] (b) 
as an example] that the K point is located at 1/6 (5/6) 
position and the K' point is located at 5/6 (1/6) po¬ 
sition on the long cutting line of ^ = 0 for metal-2p 
(metal-2m) nanotubes. Note that the positions of the 
K and K' points on the long ID BZ are the opposite to 
these on the short ID BZ. For the metal-1 nanotubes, 
the cutting lines of p = ±N/3 [p = T-N/3) pass through 
the K and K', respectively, for dx = 2 {dx = 1) where 
dx = Taod[{2n + m)/d, 3)], as shown in the previous work 
for the conventional short cutting lines 1 ^ Even though 
p = ±N/3 may exceed the range 0 < p < d — 1, the ex¬ 
pression is convenient because the K and K' points are 
mapped onto the center of the long cutting lines, which 
corresponds to ID wavenumber of fc = 0. Other choices, 
e.g. p = mod[{2n + m)/3,d] and p — mod[(2m-|-n)/3,d] 
given in Ref. [s^, may map the K and K' points away 
from the F point of ID BZ [see the longer cutting lines 
of p = 2 and p = 1 in Fig. [T] (a) for the (n, m) = (6, 3) 
metal-1 nanotube]. 


III. NUMERICAL CALCULATION 


which corresponds to shortest distance between two A 
(B) atoms in the axis direction, because of the definition 
of Note that the inversion of the range of k/2'K 

is equal to a^. Because each cutting line defined by Eqs. 
o and 0 is equal to NK 2 /d, all independent k states 


The two valleys K and K' are decoupled for the finite- 
length metal-1 nanotubes with the Cd rotational symme¬ 
try, since the two valleys belongs to states with different 
orbital angular momenta. For this case, energy levels 
show nearly fourfold degeneracy and its small lift by the 
spin-orbit interaction, as will be confirmed in numerical 
calculation in the Appendix ffiu On the other hand, the 
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two valleys can couple for the finite-length nietal-2 nan¬ 
otubes even both ends keep the Cd rotational symmetry. 
Here we perform numerical calculation of finite-length 
SWNTs to investigate the valley coupling for the metal-2 
nanotubes. Vernier-scale-like spectrum will be shown for 
an orthogonal-shaped boundary. Nearly fourfold degen¬ 
eracy and its small lift, which is not due to the spin-orbit 
interaction for shorter nanotubes, will be shown for a so- 
called minimal boundary. Vernier-scale-like spectra for 
an armchair nanotube and a capped metal-2 nanotube 
will also be shown in the Appendix TO 21 

The numerical calculation is done using the extended 
tight-binding method,— in which tt and a orbitals at each 
carbon atom is considered, and the hopping and overlap 
integrals between the orbitals are evaluated from the ab 
initio calculation^ for interatomic distances of up to 10 
bohr. Since the systems we focus on are the finite-length 
nanotubes, the electronic states are calculated by solv¬ 
ing the generalized eigenvalue problems with bases being 
from all orbitals in the systems. The optimized geometri¬ 
cal structure given by the previous energy band calcula- 
tion^i is utilized for determining the positions of carbon 
atoms. Three-dimensional structure is taken into account 
in the calculation, therefore the curvature effects^li^^ are 
automatically included. Spin degrees of freedom, and 
the atomic spin-orbit interaction Vso = 6 meV on each 
carbon atom^— are also taken into account. A tiny mag¬ 
netic field [B = 10“® T, corresponding spin splitting is 
^ 10“® meV) parallel to the nanotube axis is applied to 
separate the two degenerate states of the Kramers pairs 
in the calculation. The tube axis (z) is chosen as the 
spin quantization axis. In the following, two types of the 
boundaries are considered. The first type of the bound¬ 
ary has a geometry constructed simply cut by the plane 
orthogonal to the nanotube axis. Here we call this bound¬ 
ary orthogonal boundary. This boundary generally con¬ 
tains Klein-type terminations at which terminated site 
neighbors two empty sites^ [see red site in Fig. [5] (a) for 
(n,m) = (7,4)]. The second type of the boundary has a 
geometry removing the Klein-terminations from the or¬ 
thogonal boundary [see Fig. 131(a) for (n,m) = (7,4)]. 
For both types, the edge has minimum number of empty 
sites [dashed circles in Figs. [2] (a) and [3] (a)]. Further, 
the number of empty sites are the same with that of dan¬ 
gling bonds for the second type. The second type of the 
boundary is called minimal boundary— In the numerical 
calculation, to eliminate the dangling bonds, each single 
dangling bond at the ends is terminated by a hydrogen 
atom, and the Klein-type termination is represented by 
two hydrogen atoms4^ Both boundaries keep the Cd ro¬ 
tational symmetry of the SWNTs. 


A. Vernier spectrum 

Figure [21(b) shows the calculated energy levels ei near 
the charge neutral point for (7,4) nanotube of 50.17 nm 
length with the orthogonal boundary for both ends. Here 
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FIG. 2. (Color online) Boundary shape, calculated energy 
levels and eigenstates for (7,4) nanotube of 50.17 nm length 
with the ends of orthogonal boundary, (a) Unfolded tube near 
the left end. The empty sites are represented by the dashed 
circles, and the carbon atoms at the boundary are marked 
by the solid circles. The red sites represent the Klein-type 
termination. The shadowed areas repeat the structure in the 
unshadowed area, (b) Energy levels ei in —35 < I < 35. 
I is the energy level index numbered in ascending order of 
the energy and I = 0 corresponds to the HOMO level, (c) 
Level separation, ei+i — as a function of 1. The levels of 
— 1 < Z < 2 indicated by red color in (b) and (c) are slowly 
decaying modes, (d)-(f) Local density (d) I = —2, (e) I = —1 
and (f) I — 3. Blue shows that for A-sublattice, and red shows 
that for B-sublattice. (g) Zone-folded intensity plot of Fourier 
transform of wavefunction on A-sublattice for each level as a 
function of k. The energy for each level is added for each 
intensity plot. The intensities for states of spin-up-majority 
are shown. They are presented by either solid or dashed lines 
in turn in increasing the energy to show them clearly. The 
blue lines show the energy band calculated under the periodic 
boundary condition. Right figure in (g) shows the energy 
levels of the even parity (blue lines) and the odd parity (red 
lines) for Vso = 0. The arrows with e (even) and o (odd) in 
(g) and in (b) show the states exhibiting intravalley coupling 
with the same parity. In (d)-(g), the components of orbital 
and spin are summed up for each site or each wavenumber. 
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I is the energy level index numbered in ascending or¬ 
der of the energy, and I = 0 corresponds to the level 
of highest occupied molecular orbital (HOMO). To show 
the degeneracy behavior, the level separation, £;+i — £;, 
which corresponds to the addition energy in the tun¬ 
neling spectroscopy measurementsis plotted in Fig. 
[H (c). The levels of — 1 < Z < 2 are slowly decaying 
modes that appear in the energy gap caused by the cur¬ 
vature of nanotube surface [the local densities of I = —2, 
I = —1 and I = —3 are shown in Figs. [2] (d)-(f) as ex¬ 
amples]. The origin of the slowly decaying modes will 
be discussed in m The level separation show oscilla¬ 
tory behavior between nearly fourfold [near Z ^ — 25 
{e ^ —190 meV) and Z ^ 13 (115 meV)] and twofold 
[Z ^ —7 (£ ^ —55 meV)] degeneracies. The behavior is 
understood by (i) the asymmetric velocities between left- 
and right-going waves in the same valley pointed out in 
the previous work,™ and (ii) the strong intervalley cou¬ 
pling. The strong intervalley coupling is confirmed by 
the intensity plot in the wavenumber as shown in Fig. [2] 
(g). In this plot, the intensities for the states of spin-up- 
majority (with spin-up polarization more than 50%, in 
the calculation the polarization exceeds 99%) are shown. 
Note that the intensity at k for the spin-up-majority and 
that at —k for the spin-down-majority in spin degener¬ 
ate levels are the same because of the time-reversal sym¬ 
metry. For this case, the intensity plots for the spin- 
up-majority and that for the spin-down-majority almost 
coincides with each other for each of degenerate levels, 
that is, the orbital state of spin-up-majority and that 
of spin-down-majority are almost the same. Each eigen¬ 
function shows the strong intensity only at the left-going 
wave at the id-valley (id'-valley) and right-going wave 
at the id'-valley (id-valley). The strong intervalley cou¬ 
pling combined with the asymmetric velocities causes the 
two types of the equal interval energy levels, 
and Ziu^V/Lnt, like the vernier scalej^ The period of 
the two- to fourfold oscillation is not constant but has 
the energy dependence, for instance, the period becomes 
longer for the positive energy region. This is because the 
velocities has energy dependence reflecting the deviation 
from the linear energy band. The velocity difference be¬ 
tween left- and right-going waves becomes smaller for the 
higher energy region in the conduction band. We can also 
see the intravalley coupling when the two levels are close 
to each other, for instance, £/ ^ 83 meV and 115 meV 
as shown by arrows in the left of Fig. [5] (g). The cou¬ 
pling occurs between the same parity states which will 
be discussed in ^IVI 



kmn 


FIG. 3. (Color online) Boundary shape, calculated energy 
levels and eigenstates for (7,4) nanotube of 50.17 nm length 
with the ends of minimal boundary, (a) Unfolded tube near 
the left end. (b) Energy levels e; in —35 < I < 35. (c) 

Level separation ei+i — ei as a function of Z. The case of 
absence of spin-orbit interaction is shown by the red cross in 
the lower panel, (d)-(f) Local density for (d) I = —4, (e) 
I = —3 and (f) I = 5. (g) The zone-folded intensity plot of 
Fourier transform of wavefunction on A-sublattice for each 
level. The intensities for the states of spin-up-majority are 
shown, (h) Enlarged plot of (g) near ei = 147 meV (upper 
panel) and ei = —142 meV (lower panel). The case of absence 
of spin-orbit interaction is shown by the red dashed lines. 


B. Valley degeneracy and lift of degeneracy 

Figure [3] (b) shows the calculated energy levels £; for 
50.17 nm length (7,4) nanotube with the minimal bound¬ 
ary for both ends. The levels of — 3 < Z < 4 indicated 
with red lines are the slowly decaying modes in the energy 
gap [the local density of Z = —3 is shown in Fig. |3] (e) as 


an example]. Above and below the energy gap, the level 
separation shows almost equal interval reflecting the lin¬ 
ear energy dispersion. Fig. [21(c) shows that each of levels 
shows nearly fourfold degeneracy, which is very different 
from that in Fig. [2] (c). As shown in the lower panel of 
Fig. 121 (c), the degeneracy is lifted of the order of meV. 
The lift of the degeneracy is clearly observed even for the 
absence of the spin-orbit interaction in some energy re- 
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gion, for instance, e ^ 100 meV. From the intensity plot 
in the wavenumber in Fig. [3] (g), each level show almost 
equal intensity between the left- and right-going waves 
in the same valley, showing the valley degeneracy. The 
enlarged plot near ei = 147 meV and ei = —142 meV in 
Fig. [3] (h) shows the effect of the spin-orbit interaction. 
The four levels near ei = 147 meV show splitting into two 
spin-degenerate levels with Ae ~ 2 meV when the spin- 
orbit interaction is absent. The splitting and the intensi¬ 
ties does not change much when the spin-orbit interaction 
is turned on. The four levels near £[ = —142 meV, on the 
other hand, show almost degeneracy when the spin-orbit 
interaction is absent. When the spin-orbit interaction 
is turned on, they splits into two spin-degenerate levels 
with the splitting Ae ~ 0.3 meV. At the same time, the 
state of spin-up-majority in lower energy shows the in¬ 
tensity only at the AT-valley {kT /27r ^ 0.3) whereas that 
in higher energy shows the intensity only at the iF'-valley 
(kT/2 tt ^ —0.3). The splitting is thus induced by both 
an intrinsic property of finite-length nanotubes with the 
minimal boundary and the spin-orbit interaction. 



FIG. 4. (Color online) Length dependence of splitting of four¬ 
fold degeneracy for (7, 4) nanotube with minimal boundary. 
Three same symbols at each length show the averaged split¬ 
ting in the energy range |£i| < 250 meV except the evanescent 
modes. The splittings are then also averaged for the three 
cases of the adjacent lengths, Lnt, Lnt + a - z , Lnt + Sa^, be¬ 
cause the splitting shows nearly three-fold oscillations with 
respect to the nanotube length in the unit of = 0.022 nm. 
Solid line shows the averaged splitting with the spin-orbit in¬ 
teraction, Fso = 6 meV, dashed line shows that without the 
spin-orbit interaction. Inset shows the splitting as a function 
of the inverse of the length. The solid square with a bar (at 
0.34 meV) shows the absolute value of the spin-splitting aver¬ 
aged in the energy range |£i| < 250 meV calculated from the 
energy band calculation. 

Fig. 0] shows the splitting of the fourfold degeneracy 
for (7,4) nanotube with minimal boundary as a func¬ 
tion of nanotube length. The splitting is averaged in 
the energy range |£i| < 250 meV except the evanescent 
modes. Further, because the splitting shows nearly three¬ 
fold oscillations with respect to the nanotube length in 
the unit of a^, the average for the three cases is also 
taken. Without the spin-orbit interaction, the averaged 
splitting decreases when the length increases with 1 / Lnt 
dependence. For the actual case of the presence of spin- 
orbit interaction, the splitting for the longer nanotubes 
200 nm) are dominated by the spin-orbit interaction, 
as shown that the splitting is asymptotic to a constant 


value (~ 0.34 meV). 


IV. EFFECTIVE ID MODEL 

The numerical calculation showed that the degeneracy 
of the energy levels, and the valley coupling in the eigen¬ 
functions, are quite sensitive to the boundary shape. To 
investigate the microscopic mechanism of the valley cou¬ 
pling for given boundary conditions, we will introduce 
an effective ID model, which extract the states near the 
two valleys remaining the atomic structure. This model 
can treat microscopic analysis of the coupling of two val¬ 
leys, which can be an advantage from the conventional 
effective mass theory which treats each valley separately. 
It will be shown that the model shows evanescent and 
traveling nature of wavefunctions, and the comparison 
of the numbers of evanescent modes and the boundary 
conditions at the ends is an important key to understand 
the behaviors of valley coupling in the standing waves 
composed from the traveling modes. Analytical form of 
the discrete wavenumbers, the length dependence of the 
degeneracy lift, and the slowly decaying modes in the 
energy gap will also be shown. 


A. Effective ID model 

Let us consider the following nearest-neighbor tight- 
binding Hamiltonian, 

3 

^ + H.C., (6) 

r 

where - is the creation operator of tt electron on A 
atom at site r, and Cg is the annihilation operator 
of TT electron on the neighbor j-th B atom (j = 1,2,3) 
at r-I- Aj, Aj is the vectors from A to nearest j-th B 
atoms (see A^ in Fig. El). The summation on r is taken 
over the finite-length SWNTs. is the hopping integral 
between A and j-th B atom which is chosen as the real 
number, and, in general, they are different from one an¬ 
other because of the curvature of nanotube surface. For 
the simplicity, we ignore the spin degrees of freedom, the 
cr-orbitals, and the hopping to next nearest neighbor and 
further sites. The valley coupling in the standing waves 
can be discussed within this simplified model as shown 
later. The asymmetric velocities'^ and the spin-orbit in¬ 
teractional are not able to be captured in this model, and 
they affect as the vernier-like spectrum and the spin-orbit 
splitting as shown in the numerical calculation in Ijllll 
The construction of an effective ID model is performed 
by projecting to an angular momentum /i, as previously 
done for the achiral nanotubesiH The effective ID Hamil¬ 
tonian for the metal-2 nanotubes is given by (see the Ap- 
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pendix [C] for the derivation), 


Substituting the following forms for the solutions, 


where 


3 

jjalbi+Aij (7) 

i j=i 


A^l — ti +t2, A£2 — —t2, A£3 — —ti, (8) 


and 


= AVcr£, (12) 

for Eqs. m and (HU, one gets the following simultane¬ 
ous equations, 

+ 72A-I*^I + 73A‘i ^ (13) 

7iA-(I‘=I-‘i) + 72AI‘=I + 73A-‘i = (14) 

V 


'y ( CA,r'^a 


X! (9) 


correspond to the /i = 0 Fourier components of the op¬ 
erators. The index i is an integer specifying the ID site. 
The summation is taken over for d A atoms of d B atoms 
which satisfy zp = azi for ae and bi, respectively. Note 
that there is a pair of A and B atoms on each £ because 
there are d pairs of A and B atoms on the same z for the 
metal-2 nanotubes. The model for (n/d, m/d) = (7,4) 
SWNTs is depicted in Fig. [5l In Fig. [5l Ai atom is 
connected to B^+e, Bf_i, Bf _5 atoms. Note that the 
Hamiltonian id^=o for (7,4) SWNT is the total Hamil¬ 
tonian since d = 1. For this case the effective model has 
the same bond connection with the original model. 



FIG. 5. (Color online) Effective one-dimensional model for 
(n/d, m/d) = (7,4) SWNTs in which (ti,t 2 ) = (5,-6). The 
solid lines show the hopping between atoms. The red lines 
denote the hopping from (and to) A atom at Asite. 


B. Modes of ID model 


Since ti-l-|t 2 | = nid + mid, there are 2{nld + mjd,) sets 
of solutions (A, 77 ) for Fqs. (fT^ and (fT4)l because each 
equation is the (n/d -I- m/d)-th order algebraic equation. 
For the case of \r]\ < 1, we can call the mode A-like 
mode because the wavefunction is polarized at A atoms. 
On the other hand, we can call B-like mode for \ri\ > 1. 
Further, we can call evanescent mode at left (right) side 
for |A| < 1 (> 1 ). The mode with |A| = |? 7 | = 1 is called 
traveling mode. Let 7 be the average of jj, 7 = 
and Sjj be the difference from the average, djj = — 

7. Because of the curvature of nanotube surface, 6^j ^ 
0. Hereafter we restrict our situation to consider the 
low energy states, \£/"i\ -C 1, and small deviation of the 
hopping integrals from the average value, Idy^/yj <C 1. 
Hereafter of this subsection we will show main results for 
the total 2(n/d -I- m/d) modes. The detail derivation of 
the modes are given in Appendix 0 

It is shown that there are {n/d+m/d—2) A-like modes 
and {n/d + m/d — 2) B-like modes in the energy range 
considered. Further, by following Appendix B of Ref. 
[ 30 I . the A-like modes are classified into |t 2 | — I evanes¬ 
cent modes at left side and — I evanescent modes at 
right side, whereas B-like modes are classified into ti — I 
evanescent modes at left side and |t 2 | — I evanescent 
modes at right side. 

The remaining four modes are the traveling modes, 
or slowly decaying evanescent modes, depending on the 
energy. For the energy outside the energy gap induced by 
the curvature of nanotube surface|e| > £gap, 
there are the following four traveling modes. 


Eigenfunctions of Eq. dZD are expressed as linear com¬ 
binations of independent modes of the Hamiltonian. The 
modes of Eq. © are classified into traveling modes and 
evanescent modes. Coefficients of the traveling modes 
(and the evanescent modes) in each eigenfunction are de¬ 
termined by boundary conditions as will be discussed in 
the subsections below ( llIV Cl - IIVEI) . Here we will show 
the independent modes of the Hamiltonian. 

Hereafter we will consider the cases of n > m > 0 for 
the metal-2 nanotubes, and then we have |t 2 | = —12 > 0 
and ti > 0. For an eigenstate 1$) = with 

energy e of the Hamiltonian ©, where |cr£) is the 7r-state 
at cr-atom (cr = A, B) on Asite, we have the following 
equations of motion, 

7l<(>A^-l-(|t2|-ti) + 72<(>A£-|t2| + 73<('Af-l-ti = £4>Bi, (10) 

7l</>B^-(|t2|-ti) + l2(/Be+\t2\ + 73</>B^-ti = Sf/Ai- (II) 


{X,ri) = , ke = T'ko + k, (15) 

for the energies 

£ = + kl (16) 

Here ke in Eq. (HU denotes four wavenumbers, where 
t' = ±1 indicates the two valleys, and k, a wavenumber 
measured from t'/cq, becomes either positive or negative, 
where fco = 27r/3 -I- fep,. [Note that t = — r'/3 corresponds 
to the index for the AT (r = I) or K' (r = —I) points, 
where ,5 = I (,5 = —I) is introduced for the metal-2p 
(metal-2m) nanotubes.] k^i and ki are the shift of the 
Dirac point in K 2 and Ki direction, respectively, at the 
t' = 1 valley, because of the curvature of nanotube sur¬ 
face. From the previous energy band calculation with the 









extended tight-binding method?^ they are evaluated to 
be 


, ^^sin36> , ^,cos3d 

Kr = -Qz/jC ^2 » fcl = QzP ^2 ^ (17) 

where dt = |C?i|/ 7 r is the diameter of nanotube, 9 = 
arccos(2n-|-TO)/2\/n2 -I- m'^ + nm is the chiral angle, and 
the coefficients are evaluated to be /3' = 0.0436 nm and 
^ = —0.185 nm. The energy dispersion of Eq. (HU shows 
the energy gap 


£gap = l7l^l^l|- 


2(1 j- 


The phase $ in Eq. m is given by 
^{r'ko + k) = t' j +arg 


7 (fc -I- ifci) 


1£ 


(18) 


• (19) 


Inside the gap, jej < £gap, there are no traveling modes, 
but four slowly decaying evanescent modes, (e*'^ 
exist for the energies 




( 20 ) 


Note that k can be either positive or negative. We have 

\ki + k\ 


\v\ = 




( 21 ) 


For the case oi h > 0, |? 7 | = y /(fci + k) / {ki — k), then 
there are two B-like modes (I 77 I > 1) which are slowly 
decaying near the left end (k > 0), and two A-like modes 
(|ry| < 1) which are slowly decaying modes near the 
right end (k < 0). The decay length is estimated to 
be ^ l/|fci|, that could be much longer than that for the 
other 2{n/d + m/d — 2) evanescent modes having the de¬ 
cay length of the order of carbon-carbon bond length. 
The slowly decaying modes appeared in the numerical 
calculation as shown in Fig. [5] (e) and Fig. [3] (e). Note 
that the appearance of slowly decaying modes for each 
sublattice depends on the sign of fci. Within the nearest- 
neighbor tight-binding model, the value might be esti¬ 
mated to be positive or negative with a small value de¬ 
pending on the application of the hopping parameters 
If ki is estimated to be a negative value, the slowly decay¬ 
ing modes could appear by applying the Aharonov-Bohm 
(AB) flux42r.^ We found that fci is a positive value in 
the previous extended tight-binding calculation^ and 
the present numerical calculation showing the slowly de¬ 
caying modes consistent with ki > 0. 


C. Eigenfunctions of finite-length ID model 

An eigenstate of the finite-length ID model is ex¬ 
pressed by a linear combination of the modes derived 


in ffiVBI Let us extract only the traveling modes given 
in Eq. m in the eigenstate. Because of the real number 
of the hopping integrals in the Hamiltonian with ^ = 0, 
the components of the traveling modes in the eigenstates 
are written by real functions of standing waves, 


(t>Ae 


4>Be 


Cr COS 

r=±l 


Cr COS 

r=±l 


{ko + kr)£ — 

(fco kr) £ -\- 


4^(A)o kr) 
2 

4^(A)o kr) 
2 


9r , 
( 22 ) 
9r , 
(23) 


where r = -|-l(r = — 1 ) denotes the branch of right-going 
(left-going) wave at ki = ko- kr is the wavenumber of 
left-going or right-going waves measured from ki = kg. 
The cosine function for r = -|-l(r = —l)is the lin¬ 
ear combination of the right-going (left-going) wave at 
ki = ko and its time-reversal state of the left-going (right¬ 
going) wave at ki = —ko. For the case that the left and 
right-going waves have the same velocity, the relation 
kr=-i = —kr=+i holds in Eqs. (1331) and (1331) because 
the left- and right-going waves which compose (/ai and 
(j)Bi should have the same energy. For the asymmetric ve¬ 
locity case,— however, they are different from each other. 
We can choose that the coefficients to be positive, Cr > 0, 
and the phase 9r to be real numbers. The coefficients, 
the phases, and the discretized wavenumbers kr are de¬ 
termined by applying boundary conditions at the left and 
the right ends. 

It should be noted that the functions of Eqs. (1331) and 
(1331) themselves do not always describe the strong in¬ 
tervalley coupling, even though they have the functional 
form of linear combination of two valley states. For in¬ 
stance, when the two valleys are completely decoupled, 
each valley state as a linear combination of left- and right¬ 
going wave in the same valley is an eigenstate. Since the 
two valley states are degenerate, a linear combination 
of the two states, c+i = c_i in Eqs. (1331) and (1331) . is 
also an eigenstate of the Hamiltonian. When the spin- 
orbit interaction is introduced, the splitting of the dou¬ 
bly degenerate energy band occurs because of the lack of 
the inversion symmetry for the chiral nanotubes— >21 The 
valley degeneracy is lifted by the spin-orbit interaction, 
and the AT-valley state with spin-up (spin-down) function 
and the IL'-valley state with spin-down (spin-up) func¬ 
tion are the set of eigenfunctions for each Kramers pairs, 
as shown in the negative energy region in Figs. [3] (g) and 
(h). For this case, Eqs. (1331) and (133|) do not describe the 
orbital states for each spin state. However, for cases that 
the spin-orbit interaction is irrelevant such as Fig. [3] and 
£/ ^ 100 meV region in Fig. |31 Eqs. (1331) and (1331) could 
express the approximated orbital parts of the eigenstates 
in the finite-length metal -2 nanotubes. 
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D. Parity symmetry 

In order to discuss the boundary conditions, let us con¬ 
sider the parity symmetry of the ID model, that is, the 
Hamiltonian is invariance by exchanging ai -O- 
where Ns is the site index at the right end. This cor¬ 
responds to an inversion symmetry of ID lattice of Fig. 
0 Using the parity symmetry, it is enough to consider 
only one of the two ends instead of both ends for the 
boundary conditions. We have the following relation on 
the wavefunctions. 


(pAi = p(j)BNs+i-e for any i, (24) 


where p = ±1 are the parity eigenvalues. Applying Eq. 
(IM)) to the standing waves in Eqs. (1^ and (1^ . it is 
shown that one of the following three conditions should 
be satisfied. 


/ c+i = 0, 

\ exp (i [(fco -b fc_i) {Ns + 1) + 26i_i]) = p, 


or 


/ c-i = 0, 

\ exp (i [{ko + k+i) {Ns + 1) + 20 + 1 ]) = p, 


or 


f Cr 

\ exp (i [(fco + fcr) {Ns + 1) + 20r]) = p, 
for both r = ±1. 


(25) 


(26) 


(27) 


The selection from Eqs. (ES]) - (HZl) depends on the 
boundary type, as is discussed below. Note that the 
parity symmetry of the ID model corresponds to the 
C 2 rotational symmetry around the C 2 axis at center 
of a carbon-carbon bond in the direction perpendicular 
to the nanotube axis in the original SWNTj^ This is 
not an inversion symmetry of the original SWNT. Unlike 
the inversion operation, the C 2 rotation change the spin- 
direction. The parity symmetry for the eigenfunctions is 
broken in the presence of the spin-orbit interaction since 
the orbital state and the spin state are coupled. This is 
the reason why the parity states for the absence of the 
spin-orbit interaction are shown in Eig. (g). 


E. Boundary conditions 

There is a variety on the geometry of the boundary 
shape. Here we explicitly consider two types of bound¬ 
ary, one gives strong valley coupling and another gives 
decoupling of two valleys. Eor a moment we exclude 
the armchair nanotubes which are classified into metal-2 
nanotubes. The first type of the boundary is depicted 
in Fig. [6] (a). At the left end, both A and B atoms 
are terminated at the same z coordinate (at £ = 1). The 
end is constructed by cutting nanotube orthogonal to the 
axis direction. This is the orthogonal boundary. Each A 


(a) 



-m/d+l 1 


FIG. 6. (Color online) Two examples of the left end classi¬ 
fied into (a) orthogonal boundary condition in which |t 2 | A 
sites and ti B sites are empty, and (b) minimal boundary 
condition in which njd K sites and mjd'Q sites are empty, 
for {njd,mld) = (7,4) nanotubes. The dashed circles repre¬ 
sent the empty atomic sites and the dashed lines represent the 
missing bonds. The solid circles represent the carbon atoms 
at the boundary, and the red circle represents the Klein-type 
termination. 


atom in 1 < £ < [£ 2 ! — fi connects to a B atom, corre¬ 
sponds to the Klein-type termination.— Each A atom in 
|£ 2 | — £1 + 1 < £ < U connects to two B atoms, and each 
B atom in I < £ < |£ 2 | connects to two A atoms. Eor this 
case, the following boundary conditions are imposed, 

= 0 at £ = -|£ 2 | + 1,-'- ,0. (28) 

= ^ at £ = —£i-I-1, • • • , 0. (29) 

The number of boundary conditions for cj)Ai is |£ 2 |, and 
that for (j)Bi is £i. As shown in llIVBl in the low en¬ 
ergy region, there are totally |£ 2 | + I relevant modes (two 
traveling and [£ 2 ! — I evanescent modes) for (pAi, and 
£1 -h 1 relevant modes for (f>Bi at the left end. Therefore, 
the wavefunctions for A- and B-sublattices can be deter¬ 
mined with an arbitrarity of the amplitude. In addition, 
as shown in the Appendix the wavefunctions of A- and 
B-sublattices share a common coefficient for this case. 
Therefore, we have two solutions for the wavefunctions 
of Eqs. (1^ and (1^ : c+i = 0 and 0_i is determined, 
or, c_i = 0 and 0+i is determined. The coefficient of 
c+i = 0 or c_i = 0 reflects only the intervalley scatter¬ 
ing at the ends. For this case, either Eq. (1^51) or Eq. 
(1^51) should be satisfied. Therefore, we get the following 
discretization for the wavenumbers, 

kr = ^ + ^kr, (30) 

where Ip is an even (odd) integer for p=l(p = —1), 
6kr = {2LTT—20r) /(A^s+1) —^0 is a small offset, an integer 
L may be chosen as L = \^9r + ko{Ns + I)]/27r, where [x] 
is Gauss’s symbol representing the greatest integer that 
is less than or equal to x. In general, 5kr=+i Skr=-i, 
therefore the energy levels are not degenerate between 
r = -|-I and r = — 1. The corresponding expression is 
used in the previous article o^'^i^^i^^~ — for the armchair 
nanotubes. Eq. (1301) simply shows the discretization for 
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the standing waves of r = 1 and r = — 1. The discrete 
energy levels have generally twofold degeneracy reflect¬ 
ing the spin degrees of freedom. When we consider the 
asymmetric velocities)^ vernier-scale-like discrete energy 
levels are obtained as shown in Fig. [2j Note that the 
integer L, then the offset 6k shows nearly three-fold os¬ 
cillations when the nanotube length Ns changes because 
L ~ [{Ns + l)/3] for fco ~ 27r/3. Energy levels for finite- 
length of Ns and Ns + 3 are almost identical while that 
of iVg -I- 1 and Ns + 2 are generally different one an¬ 
other, which is confirm in our numerical calculation (not 
shown). Because the analysis above relies on the low en¬ 
ergy condition, deviation such as the intravalley coupling 
in the same parity state as shown in Fig. [2] could occur 
for a larger energy region. 

The second type of the boundary is depicted in Fig. [B] 
(b). Each A atom \iv\t 2 \ — ti + l < ^ <ti connects to two 
B atoms in the body, and each B atom in 1 < £ < |t 2 | 
connects to two A atoms in the body. This boundary is 
the minimal boundary in Fig. |31 The following condi¬ 
tions are imposed for the wavefunctions, 

= Q at £ =-|t 2 | + 1, • • • , 1 ^ 2 ! - ti- (31) 
(/)B^ = 0 at £ = —77i/c?-I-1, • • • , 0. (32) 

The number of conditions for (pAt is n/d, and that for (p-Qi 
is m/d^ Because the number of boundary conditions 
for (j>Ae is larger than or equal to the number of relevant 
modes of A-sublattice at the left end, n/d > \t 2 \ + 1, the 
standing wave for A-sublattice should be zero at the left 
end, which corresponds to “fixed boundary condition” for 
the standing waves. Then we get the condition c_i = c+i 
so that the envelope function of the A-sublattice becomes 
zero at the left end. This condition, c_i = c+i, reflects 
that only the intravalley scattering occurs at the ends. 
We also have that the phase difference 0+i — d_i is fixed 
to be 7r/2 in the linear dispersion region in which the 
relation $(fco — fc) = ^{ko + k) + tt holds. Therefore, 
the wavefunction of A-sublattice vanishes at the left end 
because of the envelope function sin(H). For this case, 
Eq. (l27l) should be satisfied for both r = 1 and r = — 1. 
We have the following discretized wavenumbers, 

k+i — fc-i _ fa -I- d+i — 9-1 , . 

2 ~ Ns + l " ^ ’ 

for both parity states p = ±1, where I is an integer. 
Note that the left hand side of Eq. (1551) is fc if fc = 
fc_|_i = — fc_i holds for the case of symmetric Dirac cone in 
which the energy is given by Eq. (HU). The corresponding 
expression is shown for the zigzag nanotubesj ^^'^'^ For a 
given I, the two parity states has the same wavenumber. 
Therefore the two states have the same energy. In larger 
energy region, we may have a parity dependent deviation 
for the phase difference from 7r/2 such as d+i — 0_i = 
tt/ 2 — dpk, where —dpk represents the deviation from 
tt/2 and ■dp is a coefficient of the deviation. For this case 
we have k = {I + 1/2)tt/{Ns -|- 1 -I- dp), then there are the 


parity splitting for the degenerate energy levels. 


Ae = I7I 


y/3a Sd 

~^Ns' 


(34) 


where dd = "dp^i — dp_i. The dependence of the inverse 
of the length on the energy splitting is consistent with 
the numerical calculation in Fig. |4l 

For the armchair nanotubes, the orthogonal and the 
minimal boundary conditions are identical. (Note that 
n/d = m/d = = |t 2 | = 1 for the armchair nanotubes.) 

Both for (pAt. and the number of boundary conditions 
is one and there are two relevant modes. Therefore the 
same discussion with the orthogonal boundary condition 
is available. 


F. Slowly decaying modes 


Finally we comment on the number of slowly decaying 
modes in the calculation. For a SWNT which is longer 
than the slowly decaying modes, the modes appear at 
the zero energy. As shown in (jlVBi the number of B- 
like evanescent modes, including the two slowly decaying 
modes, at the left end is -I-1 for the metal-2 nanotubes 
with = 0 states. For the orthogonal boundary, the 
number of boundary conditions for the B-sublattice at 
the left end is ti Therefore, the number of independent 
evanescent modes is ti -I-1 — = 1 for each spin and each 

end. There are total 4 independent evanescent modes 
for both ends orthogonal boundary. Since the longest 
evanescent mode dominates in each eigenstate, 4 slowly 
decaying modes appear. 

For the minimal boundary, the number of boundary 
conditions for the B-sublattice at the left end is m/d, 
then the number of the independent evanescent modes is 
ti -I- 1 — m/d = {n — m)/3d + 1 for each spin and each 
end. Since the number of slowly decaying modes is two 
in the independent evanescent modes, 8 slowly decaying 
modes appear for both ends minimal boundary, as shown 
in Fig. [51 The remaining independent evanescent modes 
with shorter decay length would also appear for the case 
of (n — to) /3d -I- 1 > 3. 

For the SWNTs shorter than the decay length of the 
slowly decaying modes, the slowly decaying modes ap¬ 
pear at finite energies to satisfy the boundary condi¬ 
tions, as shown in Figs. [5] and [5] Note that the A- 
like evanescent modes do not appear at the left end be¬ 
cause the number of A-like evanescent modes, |t 2 | — 1, 
is smaller than that of the boundary conditions for the 
above boundaries, |t 2 | or n/d. 


V. CONCLUSION 


In summary, we studied the discrete energy levels in 
the finite-length m-SWNTs. For the metal-1 nanotubes 
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with the Cd rotational symmetry, the two valleys are de¬ 
coupled in the eigenfunctions because they have differ¬ 
ent orbital angular momenta. The energy levels have 
nearly fourfold degeneracy and the spin-orbit interac¬ 
tion lifts the degeneracy in the order of meV. For the 
metal-2 nanotubes, on the other hand, the two valley 
states have the same orbital angular momentum and they 
are strongly coupled for the orthogonal boundary and 
the cap-termination as well as the armchair nanotubes. 
The energy levels showed the vernier-like spectrum re¬ 
flecting the asymmetric velocities and the strong val¬ 
ley coupling. For the metal-2 nanotubes with minimal 
boundary, nearly fourfold behavior was observed, reflect¬ 
ing nearly decoupling of two valleys. For this case, the 
parity splitting overcomes the spin-orbit splitting for the 
short nanotubes. The effective one-dimensional model 
explained the coupling of the two valleys, appearance 
of the slowly decaying modes caused by the curvature- 
induced shift of the Dirac point, the length dependence 
of the parity splitting. The spectrum types for nanotube 
types discussed in this paper are summarized in Table ID 

TABLE I. Summary of the degeneracy for the finite-length m- 
SWNTs with Cd rotational symmetry around the tube axis, 
and with C 2 rotational symmetry around the axis perpendic¬ 
ular to the tube axis. Metal-1 nanotubes satisfy the relation 
= d, and metal-2 nanotubes satisfy the relation du = 3d. 
OB means the orthogonal boundary, MB means the minimal 
boundary, and SO means spin-orbit. The armchair nanotubes 
with armchair edges are categorized in Metal-2 (OB). Chiral¬ 
ity dependence of the spin-orbit splitting for each valley is 
gjyen in Ref. M 

Type Spectrum 

Metal-1 nearly fourfold, SO splitting 

Metal-2 (OB) vernier-scale-like 

Metal-2 (MB) nearly fourfold, SO splitting, parity splitting 


In this paper we restrict the finite-length SWNTs with¬ 
out any external field. The finite-length effects studied 
in this paper will directly checked in devices made of 
all nano-carbons including SWNTs. In the many exper¬ 
imental setup with the metal-gated SWNTs, additional 
potentials could be created to confine the electrons inside 
the tubes. However, the nanotube ends may still influ¬ 
ence to the electrons in a center region because of the 
Klein tunneling through the barriers^ or for the nan¬ 
otubes on the metal electrodesi^ Furthermore we did 
not consider effects of Coulomb interaction between elec¬ 
trons, which can be the same order of the level spacing 
of the single-particle energy. One of the main feature 
could be captured within the so-called constant interac¬ 
tion modelji which simply increase the lab el-independent 
constant term in the addition energy. The degeneracy 
behavior studied in this paper could be useful for under¬ 
standing the orbital-related correlated electrons such as 
the Kondo effect^i^^— and the Pauli blockade<^ The in¬ 
teraction could also cause intervalley scattering as well 


as intravalley scattering, and might affects especially on 
the valley decoupling features. The effective ID model 
derived in this paper could be a lattice model to treat the 
Coulomb interaction for the given geometry of the chi¬ 
rality with the boundary for the recent observed ID cor¬ 
related electron effects such as the Mott insulator— and 
the Wigner crystallization— The effects of additional 
potential effects as well as the Coulomb interaction in the 
finite-length SWNTs with boundaries should be clarified 
in future study. 
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Appendix A: Long cutting lines passing through K 
and K' points for metal-1 and metal-2 nanotubes 

Here, we will show the fact that a long cutting line 
passes through both K and K' points for the metal-2 
nanotubes, whereas cutting line passes only through K 
or K' points for the metal-1 nanotubes. 

If a cutting line passes through both K and K' points, 
the cutting line should also passes through a F point 
because the K' point sits on the opposite side of K point 
with respect to the F point. Therefore, the cutting line 
should be /r = 0. If the /r = 0 cutting line does not pass 
a K point, there is no cutting line which passes both K 
and K' points for a given (n,m). 

The condition that the n = 0 cutting line passes 
through a K point is expressed by, 

(a + + jibi + 72^2, (Al) 

where Fi^ = (26i -|- h2)/3 is the vector from F point to 
K point in a hexagonal BZ, a, /3, ji and j 2 are integers. 
Here /? is introduced as follows; 

{ 0 for metal-1, 

-t-1 for metal-2p, (A2) 

— 1 for metal-2m. 

In Eq. dAH), we used the already known fact that a K 
point is mapped onto the center of a short cutting line 
for the metal-1 nanotubes, whereas it is mapped onto 
5/6 (1/6) position of a short cutting line for the metal- 
2p metal-2TO nanotubes— It should be noted that the 
following relation holds for the metal-2 nanotubes, 

mod (^, 3 ) = mod 3^ , (A3) 
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because of the relation (n — m)/3d = [(2n + m) — {2m + 
n)]/dii = — In order to satisfy Eq. (lAlL ji, j2 and 

a should satisfy the following equations. By comparing 
the coefficients of 6 i and 62 in Eq. m, we have, 


jin + j 2 m = -^{2n + m), 

(A4) 

[2 \N /3 

a = 0 +-A-T- 

Vo J m 0 

(A5) 


Eor the case of metal-2 nanotubes, it will be shown 
that the following ji and j 2 satisfy the conditions Eqs. 
(IA4[) and (IA5I) . 


1 / 

m \ 

(A6) 

= -3 1 

.^7 + ’ 

1 / 


(A7) 

= -3 1 

~^d + ■ 


Note that the right hand sides of Eqs. (|A6p and (lATl) 
are integers because of Eqs. (0 and (jA3l) for the metal-2 
nanotubes. By substituting Eq. (jA6l) for Eq. (IA5I) . a is 
given by 


Of = — 


3 



(AS) 


Because N = 2(n^ + m"^ + nm)/dfi = [(2m -I- n){2n + 
m) — 3nm\/dii and dn = 3d for the metal-2 nanotubes, 
the following relation folds. 


N nm 

=-3tif2-^. 

Using this, Eq. (IA8I) becomes 


a = 




(A9) 


(AlO) 


The right hand side of Eq. (jAlOll is an integer because 
of the relation mod(nm/cP, 3) = 1, which can be derived 
from Eqs. (O and (IA3I) . Therefore, Eq. m is satisfied 
by the set of integers of Eqs. (IA6I) , (IA7I) and (lAlOl) for the 
metal-2 nanotubes. Erom the left hand side of Eq. (ED) 
and Eq. (lASp . it is shown that the K point is located at 
1/6 (5/6) position of the longer cutting line defined by 
Eqs. (H]) and ([2]) with /r = 0 for the metal-2p (metal-2m) 
nanotubes. 

For the case of metal-1 nanotubes, dividing Eq. ED 
by d, one gets. 


Appendix B: Numerical calculation for metal-1, 
armchair and capped metal-2 nanotubes 

Here we will show some numerical calculation of the 
finite-length metal-1, armchair, and capped metal-2 nan¬ 
otubes using the extended tight-binding model. As ex¬ 
pected, the metal-1 with Cd rotational symmetry will ex¬ 
hibit nearly fourfold degeneracy and its lift by the spin- 
orbit interaction. On the other hand, the armchair and 
the capped metal-2 nanotubes will exhibit the vernier¬ 
like spectra as well as Fig. [2] 


1. Energy levels for metal-1 nanotubes 
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FIG. 7. (Color online) Boundary shape, calculated energy 
levels and eigenstates for (6,3) nanotube of 50.15 nm length 
with both ends minimal boundary as shown in (a). (b) Energy 
levels ei in —35 < Z < 35. (c) Level separation, ei+i — ei, as 
a function of 1. The dashed lines in the lower panel show the 
spin-orbit splitting for corresponding energy calculated by the 
energy band calculation, (d)-(f) Local density for (c) I = —4, 
(d) I = —3 and (e) I = 5. 


41- 


. m 1 
■ J2-7 = ■ rZ 2- 


(All) 


Because t 2 is not a multiple of 3 for the metal-1 nan¬ 
otubes, any integers of ji and j 2 cannot satisfy Eq. 
(lAllI) . Therefore, there is no cutting line passing through 
both K and K' points for the metal-1 nanotubes. 


Figure|7](b) shows the energy levels for (6,3) nanotube 
with 50.15 nm length with the minimal boundary for 
both ends keeping C 3 {d = 3) rotational symmetry. The 
left end is depicted in Fig. [7] (a). There are the slowly 
decaying modes (—3 < Z < 4) in the energy gap between 
e /=_4 = —44 meV and e ;=5 = 53 meV. Above and below 
the energy gap, the level separation shows almost equal 


















































13 


interval reflecting the quantization of the linear energy 
dispersion. Fig. [7] (c) shows that each levels show nearly 
fourfold degeneracy. The levels show complete fourfold 
degeneracy for the case of absence of spin-orbit inter¬ 
action. The spin-orbit interaction lifts the fourfold de¬ 
generacy as expected in the energy band calculationj^i 
Other metal-1 nanotubes, for instance, metallic (9,0)- 
zigzag nanotubes, also show similar behaviors with Fig. 
[3 nearly fourfold degeneracy and lift of the degeneracy 
by the spin-orbit interaction (not shown). In addition, 
finite length (9,0) nanotubes show the edge states decay¬ 
ing in the length scale of carbon-carbon bond, not on but 
below the charge neutral point, as the effect of the next 
nearest-neighbor hopping processi^ 



(b) 

J, 

















FIG. 8. (Color online) Boundary shape and calculated energy 
levels ei for (6, 3) nanotube of 50.63 nm length, (a) Unfolded 
tube near the left (lower) and right (upper) ends. The upper 
end is cut along the na\ and ma^. Note that both ends are 
classified into the minimal boundary, (b) Energy levels ei in 
—35 < I < 35. (c) Level separation, si+i — ei, as a function 
of L 

When the rotational symmetry at the end is broken, 
the valley degeneracy and its lift by the spin-orbit inter¬ 
action will not be clearly observed. Figure H] shows the 
energy level for (6,3) nanotube with 50.63 nm length, in 
which the one side of the end is cut along the nai and 
77102 and the other side is the same with Fig. [3 The 
system loses the C 3 rotational symmetry which is pos¬ 
sessed in the case of Fig. [3 The large lift of the fourfold 
degeneracy, for instance, Ae > 6 meV in e ^ 200 meV, is 
not caused by the spin-orbit interaction, but due to the 
mixing of the two valley degrees of freedom by breaking 


of the rotational symmetry. 

In the end of this subsection, we comment on val¬ 
ley mixing effect by the spin-orbit interaction. Strictly 
speaking, the spin-orbit interaction could mix the states 
in two valleys because spin-up states in /r = 1 cutting 
line, which passes the K' point, and spin-down states in 
/i = 2 cutting line, which passes the K point, have the 
same total angular momentum, 3/2 (see Fig. [3 (a) for 
the cutting lines). Mixing of the two valleys may give 
an additional effect from the spin-orbit splitting, such as 
the parity splitting in Fig. [3] or the splitting in Fig. [5] 
for the absence of C 3 rotational symmetry. Such valley 
mixing effects, however, seem to be irrelevant in Fig. [3 
which simply shows the splitting as expected from the 
band calculation. 


2. Vernier spectra for (6,6)-armchair and capped 
(7,4) nanotubes 



(b) 

1 1 . 


















-0.35 -0.3 0.3 0.35 


kT/2j[ 

FIG. 9. (Color online) Boundary shape, calculated energy 
levels and intensity plot in wavenumber for (6, 6) armchair 
nanotube of 50.05 nm length, (a) Unfolded tube near the 
left end. (b) Energy levels in ei =_35 < e < £ 1 = 35 . (b) Level 
separation, £i+i — ei, as a function of 1. (c) Intensity plot 

of Fourier transform of wavefunction on A-sublattice for each 
level. The intensities for the states of spin-up-majority are 
shown. The blue lines show the energy band calculation under 
the periodic boundary condition. 

Figure [3 shows the energy levels for (6,6)-armchair 
nanotube with 50.05 nm length. The energy levels show 
the similar behaviors with the vernier-like spectrum in 










































































































































14 


Fig. This can be understood by (i) the strong inter¬ 
valley coupling and (ii) the asymmetric velocities, as well 
as Fig. [2] For the armchair nanotubes, it can be shown 
using the effective ID model that even channel from 
{ae + hi)l\f2 states and odd channel from [ai — &^)/-\/2 
states are decoupled each other for this boundary, and 
the even (odd) channel has the energy band with left¬ 
going states at the iF-valley (iF'-valley) and right-going 
states at the iF'-valley (iF-valley). Therefore, no intraval¬ 
ley mixing is seen in Fig. (c). The period of the two- 
to fourfold oscillation is not constant but has the energy 
dependence, for instance, the period becomes longer for 
the positive energy region. This is because the velocities 
has energy dependence reflecting the deviation from the 
linear energy band. The velocity difference between left- 
and right-going waves becomes smaller for the higher en¬ 
ergy region in the conduction band as well as the (7,4) 
nanotube. 




kmit 


trum for a capped nanotube, which would be more abun¬ 
dant than the orthogonal boundary containing the Klein- 
terminations. Figure [13 shows the calculated energy 
levels for both-side-capped (7,4) nanotubes of 49.56 nm 
length. In the calculation, the cap structure is formed 
by using the graph theory ;^^'^^ in which the cap region 
is defined outside the end cut along the nai and ma 2 . 
There are two cap obeying the isolated pentagon rule for 
(7,4) nanotube, one has 55 and the other has 57 car¬ 
bon atoms at the cap region. The cap with 57 carbon 
atoms is used in the calculation. To connect the cap 
and the body smoothly, the structure is optimized by 
the molecular mechanics method with the Universal force 
field (UFF) in the Gaussian program<^ Even the opti¬ 
mization method gives less accuracy on the electronic fine 
structure such as the curvature-induced energy gap and 
the spin-orbit interaction, we could discuss the intervalley 
coupling in the eigenfunctions from the semi-quantitative 
calculation. The vernier-like spectrum is seen in the cal¬ 
culated energy levels. As shown in the plot of Fourier 
transform, in general, each level is formed from a left¬ 
going wave of one valley and a right-going wave of an¬ 
other valley, that corresponds to \cr/c-r \ <C 1 (r = 1 or 
— 1) in Eqs. (1221) and (1231) reflecting the strong interval- 
ley coupling. For closer two levels, the intravalley mixing 
between the same parity states is also seen, for instance, 
£i ^ —208 meV, 181 and 213 meV. The vernier-like spec¬ 
trum similar with Fig. [10] is also obtained for another 
cap obeying the isolated pentagon rule with 55 carbon 
atoms in the numerical calculation (not shown). 

The strong valley coupling for the case of cap- 
termination may be understood in the ID model as fol¬ 
lows. The boundary conditions at the left end can be 
given by (j)Ai = 4>ct^ for n sets of (^a, 4) and (j)Bi = (t>cic 
for m sets of {£a,^c), where Ia (^b) is the index of the 
coordinate of A-sublattice (B-sublattice) and ic is that 
of cap. For the amplitudes 4>cict there are Nc equations 
of motion. In general, (f’dc shows oscillation in the length 
scale of carbon-carbon bond. Therefore, one also expects 
the fast oscillations for the wavefunctions of A- and B- 
sublattices, which can be formed under the strong inter¬ 
valley coupling of \cr/c-r\ 1 for r = 1 or —1. 


FIG. 10. (Color online) Boundary shape, calculated energy 
levels and intensity plot in wavenumber for (7,4) nanotube of 
49.56 nm length with both-ends-capped. (a) 3D and (b) 2D 
representation of the cap structure. In (b), the solid circles 
indicate the A and B atoms connected to the cap region, and 
the open circle indicate the carbon atoms in the cap region 
at the left end. (c) Energy levels ei in —35 < I < 35. (d) 
Level separation ei+i — ei as a function of 1. (e) Intensity 

plot of Fourier transform of wavefunction on A-sublattice for 
each level. Right figure in (e) shows the energy levels of the 
even parity (blue lines) and the odd parity (red lines) for 
Vso = 0. The arrows in (c) and (e) show the states exhibiting 
intravalley coupling. 


We show another case exhibiting the vernier-like spec- 


Appendix C: Derivation of effective ID Model 


Here we will show the detailed of the derivation of the 
effective ID model given in lUV AI 

In the cylindrical coordinate system in which the tube 
axis coincides the z axis, the components of Aj appearing 
in Eq. ([S]) (A^ is depicted in Fig. |TT]) is given by 


A9l = TT 

A02 = ~7r 

A03 = —TT 


n + m 


rfi + m? + nm ’ 
m 

-|- -|- nm ’ 

n 

-|- m^ nm ’ 


(Cl) 

(C2) 

(C3) 
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and 


Azi = —a 
Az 2 = a 


n — m 
' M ’ 

2 n + m 


Az 3 = -a. 


3d ’ 
2 m + n 
3d ■ 


(C4) 

(C5) 

(C 6 ) 


Let us consider the following Fourier transform in the 
circumference direction for the operators, 








(C7) 


where ct = A, B, and z = azi denotes the lattice position 
of A atoms on the nanotube axis, ^ is the integer for the 
lattice position. The summation on r is taken place for 
a given z. The inverse Fourier transform is given by 




1 

7^ 


d-l 


/i—0 


(C 8 ) 


Substituting Eq. (IC 8 I) for the Hamiltonian Eq. ([B]), the 
Hamiltonian can be decomposed into projected Hamilto¬ 
nian for /i-th angular momentum 


H = Y,H^, (C9) 

3 

+ H.C., (CIO) 

z j=l 


By selecting the /i values in which H^j contains states near 
the Fermi energy, one gets an effective ID Hamiltonian. 


ed,l2 


FIG. 11. Coordinates for (7,4) nanotube. The dashed lines 
show the interval of the ID lattice constant = Td/N. For 
this case, d = 1, T = v^a and A = 62 . 

Because of our main interest, let us focus on the metal- 
2 nanotubes. For the metal-2 nanotubes, Azj can be 
expressed by 


Azi = Qzih -1- £ 2 ), 

(Cll) 

Az2 = -azt2, 

(C12) 

Az3 = -ttzh. 

(C13) 


It should be mentioned that there is a B atom at pCh /da 
from each A atom, that is, there are d pairs of A and 
B atoms on the same z for the metal-2 nanotubes. [For 
example, see Fig. [TT]for the coordinates of (n, m) = (7,4) 


metal-2p nanotube {d = 1).] We consider only i7^=o and 
use the simplified notation, 

CA(/i=0,z=a^^) 0^1 CB(iJ,=0,z=a^£) b^. (C14) 

and A£j = Azjjuz, then we get Eq. ([7]) for the effective 
ID Hamiltonian for the metal-2 nanotubes. 


Appendix D: Mode analysis of effective ID model 

Here we will discuss the detail derivation of the modes 
of Eqs. (fT3|l and (fT4)l . 

As a general property, when (A, rj) is a set of solution of 
Eqs. (ESI) and (ITi)) . (1/A, l/ry) is another set of solution 
since the equations are equivalent by changing (A, 77 ) O 
(1/A, l/rj). For a solution (A, 77 ), the complex conjugates, 
(A*, 77 *) is also another set of solution since and e are 
real numbers. 

We will start from the traveling modes (|A| = I 77 I = 1) 
in the solutions of Eqs. (H and O- For the traveling 
modes, Eqs. m and o are equivalent as mutual com¬ 
plex conjugate. Let us first consider the flat graphene 
case, 6 jj = 0. By noting the fact that the left hand side 
of Eq. (fT51) is zero for A = because of the relations 

Eqs. m and (IA3I) . we can expand at the wavenumbers 
ki = ±27r/3 which correspond to the two Dirac points. It 
is shown that there are four traveling modes for a given 
energy; (A, 77 ) = where ke = r'27r/3-|-fc' de¬ 

notes four wavenumbers for a given energy. We restricted 
in the linear dispersion regime, in which the energy and 
k' has the relation, 


e = ±| 7 ||E|A:'|. (Dl) 

The phase in Eq. m is given by 

* ("'f + ‘') =(f +'®'') + "* (7) ’ 

/I = ±1 is the value defined in Eq. (IA2|) for the metal- 
2 nanotubes, and 6 is the chiral angle. Because of the 
time-reversal symmetry, we have $(fc^) = — $(—fc^). 

When Sjj ^ 0, the above analysis should be modified 
as follows. For a state with wavenumber ki = t'2'k/ 3 + k', 
by considering the contribution of lowest order of fc' and 
/y, Eq. ((T5)) is written as, 

— iye"^ i(^t2-i-/36') 7^ 

2 az 

Here k^ and ki satisfy the following relations, 
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/cr and ki relate to the shift of the Dirac point in iiTi 
and K 2 direction, respectively, from the K or K' points 
in 2D BZ because of the curvature of nanotube surface. 
By comparing with the energy band calculation with the 
extended tight-binding method;^ they have the relations 
in Eq. dm)- From Eq. (|D3I) . it is shown that there are 
the following four traveling modes, 

(A, ry) = , ke = T'ko + k, (D 6 ) 

for the energy outside the energy gap, |e| > £gap, where 
the energy gap 


Egap = l7l|^l^l| (D7) 

is induced by the curvature of nanotube surfaceJ^* 2 iii 2 
After Eq. (ID 6 p . k is the wavenumber measured from 
r'fco = t' (27r/3 -I- A:r), which is the bottom (top) position 
of the conduction (valence) band. The energy and k has 
the following relation 




(D 8 ) 


and the phase is given by 
^{r'ko+k) =t' + I3e\+aig 


'y{k + ifci) 


le 


• (D9) 


Inside the gap, jej < £gap, there are no traveling modes, 
but four slowly decaying evanescent modes exist. Near 
the band edges, the modes can be analyzed by changing 
{k' — t'/cr) —?> in in Eq. (ID3I) and in a pair equation 
derived from (na. Then it is shown straightforwardly 
that the modes (e*"^ ^““'', 77 ) have the energy 

Note that k can be either positive or negative. From Eqs. 
(ID3I) and (IDIOI) . we have 


\ki 




(Dll) 


The remaining 2{nld+mld) — A modes can be classified 
into {njd + mjd— 2 ) A-like modes and {n/d + m/d — 2) 
B-like modes. For the A-like modes, using je/yl <C 1, 
l^7i/7l ^ 1 9 ,nd \r]\ < 1, Eqs. (fT^ and (03 can be 
simplified as 

X^/d+m/d ^n/d + 1 = 0, (D12) 

£ 1 
7A-(l*2l-*i)-hAl*2l-HA-‘i' 

The first equation (ID12I) has the roots A = 

and the properties for these modes have already been 

captured as the traveling modes or the slowly decaying 


modes in the previous paragraphs. For the remaining 
{n/d+m/d—2) solutions, the analysis given by Akhmerov 
and Beenakker— is applicable. Note that Eq. (ID12P for 
A' = xn/d+m/d ig Eq^3.7a) in Ref. IM. 

By following Appendix B of Ref. it is shown that 
there are \t 2 \ — 1 roots which satisfy |A| <1 (evanes¬ 
cent modes at the left side), and — 1 roots which 
satisfy |A| >1 (evanescent modes at the right side) for 
Eq. (ID12I) . The second equation (ID13I) determines ry for 
each A. The similar discussion can be done for the B-like 
modes with the equations which are given by changing 
n o m (then ti o \t 2 \) and ry —^ l/yy in Eqs. (ID12I) 
and (ID13I) . Then it is shown that there are — 1 modes 
which satisfy |A| < 1, and there are |t 2 | — 1 modes which 
satisfy |A| > 1 for the B-like modes. 


Appendix E: Traveling modes for A- and 
B-sublattices for a given boundary 


Here we will show that one of the following relations 
should be hold for the orthogonal boundary: c+i = 0 
and 0-1 is determined, or, c_i = 0 and 9+i determined 
in Eqs. (1^ and (^3 . 

Let us first show a useful relation. The left hand side 
of Eq. 113 is rewritten as 


t2-l 




AI*2| 


n n n 


t'=± 


mi—1 
* 2-1 


m2 = l 


=T n (^-v) n 


r'—:t 

t 2 -l 


mi — 1 


1 1 \ 

n T-i^) n (A-A>.), 


mi — 1 




(El) 


where Xr' = , X ^^, and X^^ are the roots of the 

left hand side, and we have |A<J < 1 , > 1 - In 

Eq. dHI), it is shown that the value rimiil 
is a real number because there is another root of A)))* 
for a complex root X/Xj^- For a traveling mode of A 
which has almost the same wavenumber with ko or — fco, 
(A — At') is a small real number. Let us ex¬ 
plicitly consider a traveling mode close to A+, (A+,e‘'®’). 
From Eqs. USD and (fETI) . we have the following relation, 
which will be used later, 


* 2-1 

n 

mi—l 


1 



(E2) 


where i? is a finite real number. 

For the limit of |e/ 7 | <C 1, the equations of motion 
for A- and B-sublattice are decoupled [see Eqs. m 
and (ITT]) ]. For this case, the wavefunctions of A- and B- 
sublattices under a boundary are determined separately, 
in general. We will show that, for the orthogonal bound¬ 
ary, the wavefunctions of A- and B-sublattices share a 
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common coefficient for the traveling modes. The wave- 
functions near the left end are written as the linear com¬ 
bination of the relevant modes as follows: 

1*21-1 

</>A^ = ^ ^ (E3) 

t'—± mi—1 

tl — 1 ^ 

4>Be = ^ CBT'e"^-I- ^ CBmaT^- (E4) 

r'=± m 2 = l ^"*2 

Note that there are two traveling modes for each sublat¬ 
tice in the limit of \e/^\ <C 1. The factors are 


attached for the convenience of the discussion. Here we 
used that each evanescent mode of B-sublattice at the left 
end has a pair with that of A-sublattice at the right end, 
that is, for each pair of roots, the relation Aa = 1/Ab 
holds where Aa is a root of Eq. 0 and Ab is the cor¬ 
responding root of Eq. 0. The coefficients are deter¬ 
mined by employing boundary conditions. If cat' = cbt' 
is satisfied for both r' = ±1, one of the following rela¬ 
tions should be hold in Eqs. (l22l) and (|23l) : c+i = 0 and 
9-1 is determined, or, c_i = 0 and is determined. 

Let us consider the orthogonal boundary expressed in 
Eqs. (1^ and (1^ . If cat' = cbt' is satisfied for both 
t' = ±l for this boundary, the following determinant 
should be zero. 


D = 


e 2 


i± 

e 2 


e^tAl^ 


g-i|A-l‘2l+i eitA:l‘=l+' 


(Af 


(A<)-I*2Hi 




e 2 


‘ = at 


e 2 




0 


-1*21 + 1 


fe)' 


e‘2 A, 


e ' 2 A_ 


(^) 


— tl + 1 




-1 


— tl + 1 


By cofactor expansion and use the relation on the Vandermonde matrix, one gets the following relation. 


D={-1) 


ti(ti-l) , |t2l(|t2|-l) 
O ' 


n (^-A^) n k.-a;;.) 


f / I 1 \ *2-1 / 1 \ i*ii-i 'I 

n ) n (A--As.)-e-» n (^-—) n (a+-a>.) 

[ mi = l ^ ^2 = 1 mi = l ^ ^2 = 1 J 


(E5) 


(E6) 


By using Eq. (IE2I) in the second line of Eq. (IE6I) . one gets should be held for the orthogonal boundary; c+i = 0 and 

D — 0. Therefore, one of the the following two relations 9-i is determined, or, c_i = 0 and 9+i is determined. 
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